NORDITA-02-42 AP 
ULB-TH/02-22 



Photon- Axion mixing in an inhomogeneous universe 



Mattias Christensson* 
Nordita 

Blegdamsvej 17, DK-2100 Copenhagen, Denmark 

Malcolm Fairbairn^ 

Service de Physique Theorique, CP225 
Universite Libre de Bruxelles, B-1050 Brussels, Belgium 

April 2003 



Abstract 

We consider the dimming of photons from high redshift type la supernovae through 
mixing with a pseudoscalar axion field in the intergalactic medium. We model the elec- 
tron density using a log-normal probability distribution and assume frozen in magnetic 
fields. Assuming the magnetic fields are produced in the early universe we are unable 
to obtain sufficient dimming in order to explain the apparent acceleration without vi- 
olating the bounds on the frequency dependence of the dimming. We also show that 
any axion mixing leading to a reduction in optical luminosities would also lead to a 
significant reduction in the polarisation of UV light from intermediate redshift objects 
which may be detected in the future. 



1 Introduction 



At the time of writing, the present understanding of the energy content of the universe 
appears to be converging towards a standard model. Observations of the cosmic microwave 
background radiation seem to suggest that the total amount of energy in the universe is such 
that there is zero spatial curvature on the largest observable scales jl]. At the same time, 
observations of the distribution of matter at smaller scales give us information about galaxy 
clustering and the time of matter-radiation equality, and tell us that the fraction of this total 
energy which is in the form of matter is only about 30% f?. This is much larger than the 

*E-mail: mattias@nordita.dk 
^E-mail: mfairbai@ulb.ac.be 



1 



fraction of the total energy which can be in the form of baryonic matter [3], and together 
with observations of the rotation curves of galaxies and velocity dispersions of galaxies within 
clusters, has lead to the concept of dark matter. However, there still remains the question 
of the remaining unaccounted for ~ 70% of the total energy density. 

Although the data set is effectively still quite small, observations of type la supernovae 
seem to have cast some light on this missing energy |^ |3] . The apparent dimness of these 
supernovae at high redshift suggests that there is a non-zero cosmological constant which 
makes up the missing energy or at least something which acts like one. If we are to interpret 
this cosmological constant as some energy field then we find that its energy density p ~ 
{lQ~^eVY is far away from any of the energy scales (i.e. the Planck scale or the scale of 
supersymmetry breaking) one would naively expect such a field to have if it were non-zero. 
It is also not clear how such a small energy scale could be stable to quantum corrections. 

Recently it has been suggested that the dimness of high redshift supernovae is not due 
to the accelerated expansion of the universe, but rather due to mixing between the photons 
coming from these objects and a pseudo-scalar axion field jB]. The Lagrangian density of 
the photon-axion system is given by 



where F^^ is the electromagnetic kinetic term and F^^y its dual, a the axion field, the 
axion mass and Ma is its coupling to the photon field. The coupling of the photon and axion 
fields in this way means that a photon has a finite probability of mixing with its opposite 
polarisation and with the axion in the presence of an external magnetic field jTj, such as 
the intergalactic magnetic field. If the magnetic field is distributed randomly along the path 
of the light this will lead to a gradual reduction in luminosity until only 2/3 of the initial 
energy of the flux remains in the form of photons. 

If the axion mass and coupling scales are chosen to be ma ~ 10~^^eV and Ma ~ 4 x 
lO^^GeV respectively jHj, the fractional amount of dimming will not vary much over the 
optical range where the supernovae are observed. This might enable one to consider models 
where the extra energy component of the universe is not a cosmological constant but a fluid 
with a different equation of state, e.g. a string network [S]. The first paper suggesting 
such a mixing broke the line of sight up into cells within each of which the magnetic field 
was oriented in a random direction |B] . Later papers included the effect of a redshift varying 
magnetic field and the presence of intergalactic electrons P[TU| (see also [H]). If the electron 
density is assumed to be approximately uniformly distributed in the intergalactic medium 
with the value fie = 1.8 ■ lO^^cm"^ the oscillations will strongly depend on the energy of 
the photons leading to possible conflict with supernovae spectral observations However, 
as pointed out in the intergalactic medium is not uniform and many regions will have 
a electron density much below the mean density. This, according to ^21; will render the 
energy dependence unobservable at present. 

In this paper we study the effects of the variation of the free electron density in the plasma 
on the photon-axion mixing. Due to the large electrical conductivity of the intergalactic 
medium the magnetic field can be considered frozen into the plasma. The magnitude of the 
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magnetic field will therefore depend on the electron density as -B oc nj . In regions with 
a low electron density the mixing probability will be suppressed because in these regions 
the magnetic field is reduced. Therefore, when both electron density and magnetic fields 
are allowed to fluctuate, the total mixing probability for a path along the line of sight is 
sensitive to regions where fluctuations have produced an electron density high above the 
mean density. It is therefore important to take fluctuations in both the electron density and 
the magnetic field into account when evaluating the total mixing probability. 

Throughout this paper we have assumed that the intergalactic magnetic fields were cre- 
ated in the early universe, i.e. at much higher redshifts than any of the objects we will be 
considering. There are some models where the intergalactic magnetic fields are produced 
at much later epochs. For example it is possible that super massive black holes may be 
responsible for the amplification of magnetic fields in the intergalactic medium at late times 
[TSj . Different results would be obtained if one were to consider such models. 



In this section we will describe the equations and assumptions we have used in simulating the 
passage of photons through the intergalactic medium. We define the background cosmology, 
describe the photon-axion mixing and then discuss the probability distribution used to model 
the electron density and magnetic fields. 

2.1 String network cosmology 

We use the standard notation for each component of stress energy in the universe such that 
the density of component x is given by Vtx = Px/Pcrit where pcrit = SHq/SttG. The equation 
of state of the component x is denoted uj^ {Px = oj^Px)- The model in [H] assumes that 
the overall density of matter VLm ~ 0.3 as suggested by observation, but that the remaining 
dark energy takes the form of a string network with a different equation of state to that of 
a cosmological constant jH]. The energy budget is therefore divided up as 



where the energy density of the string network is denoted Vtg. The energy density of the 
string network allows the universe to possess a flat geometry whilst the equation of state of 
the string component means the Friedman acceleration equation is identical to that of an 
open matter dominated universe. 

The distance / travelled by a photon reaching us today emitted at a redshift within 
this cosmology is 



2 Simulations 
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Supernova results suggest an ^Im = 0.3, f^A = 0.7 universe and one can evaluate the cor- 
responding equations for such a cosmology. In order for the luminosity distance of such 
a universe to be mimicked by dimming in a string network universe as described in equa- 
tion (j2I), one would require a dimming at redshifts of 2; ~ 0.5 and 1 of about 12% or 15% 
respectively. 

In order to see if it is possible to obtain such a dimming we need to investigate the 
photon-axion mixing. 

2.2 Photon-axion mixing 

We need to investigate how an ensemble of photons of different polarisations will evolve in 
the presence of external electro-magnetic and axion fields. In order to see fully how the 
flux and polarisation coming from a source will change we will need to start with an initial 
density matrix which encodes the information about the polarisation of the source. However, 
we will introduce the mixing equations using normal vector potentials. 

Our analysis of the mixing follows closely the analysis of [0] and ^U]. Like them we ignore 
the effects of vacuum polarisation J3] and the birefringence of the intergalactic medium in 
magnetic fields. 

We write A± and A\\ for the polarisations of the photon perpendicular and parallel to the 
magnetic field respectively. Then since the refractive index is close to unity we can linearise 
the equations of motion for the axion a and the two polarisations of the photon 

= (4) 

where u is the energy of the photon and the mixing matrix Ad is given by 

/ Ap \ 
M=\ Ap Am \. (5) 
\ Am A^ / 

Here the quantities Am, Ap and A^ are defined as 

— -(lo^)X^) 

-3.6 X 10-- ('-^) (^^). 

\ 00 J V10-7cm-3y 

where B± is the magnetic field strength perpendicular to the direction of propagation of the 
photon. 
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The eigenvalues of the mixing matrix (jSJ are given by 

(Ap + A, 



Ai = A, 



and 



A 



2,3 



A^ a2 
4 4 



1/2 



(6) 



where Ai is trivial since the photon can only mix with the polarisation state of the photon 
parallel to the magnetic field. The mixing angle which diagonalises the bottom right hand 
four components of (0) is 

^" (7) 



tan26' 



A„ - A, 



We name the rotation associated with this mixing angle Ui 

Ui = 



1 
COS 6 sin ( 
— sin 6' cos ( 



so that we can write the evolution equation for the three states in the form 

= uIdUi 





Here the matrix D is the diagonal matrix consisting of the propagating eigenstates 



D = diaKfe-^("+^^)*,e'^("+^^)\e-'(^+^3)'l 



(9) 



(10) 



and equation © gives the evolution of the A\\ and a amplitudes over a single domain. 

In our analysis, the photon/axion states pass through many regions with randomly ori- 
entated magnetic fields. Because of this, Ay in one region or cell is not the same as A\\ in 
the next and we have to perform an additional rotation within each cell to take this into 
account. It is convenient to express the evolution directly in terms of the A^, Ay where x 
and y are fixed directions perpendicular to the direction of propagation. 
We use the name U2 for the rotation between the states A^j-^y and Aw 
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where ip is the angle between the A\\ and Ay components. We can then write 

= UlulDUiU2 

and we introduce the notation for this matrix 

T = ululDUiU2. 





(12) 
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Now the unitary matrix T is a function of ng, B± and ip so it is different for each cell. We 
therefore call the T corresponding to the nth cell T„, the time of flight after n cells t„ and 
the photon or axion states after n cells A{tn) and a(t„) respectively. We can then write the 
evolution after n domains as 



a{tn) 




TjiTn-l-.-Ti 




(14) 



This shows that in order to work out the effect of mixing in each of the domains along the 
path of the photon/axion states we can simply multiply the matrices Tn from each domain 
along that path. We obtain the total amount of axions and photons after a time t„ by 
calculating the total conversion probability at the end of the path. 

Strictly speaking, this approach is only valid for a pure initial state. For unpolarised 
states where the initial state vector has contributions from two perpendicular polarisations 
one should evolve the full density matrix keeping track of the off-diagonal elements. How- 
ever, since the magnetic field direction and strength is randomly generated in each cell, the 
information in the off-diagonal parts of the density matrix will not be important in the 
determination of the expected amount of flux and polarisation after a certain time. 

We have verified this in our simulations by tracking the evolution of the full density 
matrix p using 

^^,p=^[M,p]. (15) 

We can use the same matrices introduced earlier 

p(t„) = T„ p(t„_i) Tt (16) 

so that the density matrix p final which gives rise to the probabilities where the wavefunction 
is broken down, i.e. where the photons are actually observed by telescopes on earth, is given 
in terms of the emitted state at the source PiniUai by 



P final 



Ttotal Pinitial ^foiai (l'^) 



where Ttotai = T„T„_i...Ti. 

For a given magnetic field strength, the mixing probability Pceii increases with decreasing 
density. Also, the energy dependence of the dimming becomes smaller with decreasing 
density. So the total amount of dimming along each line of sight depends critically upon the 
behaviour of both the magnetic field and the electron density along each line of sight. This 
is what we turn to next. 



2.3 Distribution of electrons and magnetic field 

Much of our treatment of the distribution of electrons will follow the investigation into 
cosmological magnetic fields of Blasi et al jT3]. The high electrical conductivity of the 
intergalactic medium allows one to assume that the number density of electrons traces out 
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the cosmic matter distribution. Therefore we need to estimate the spectrum of density 
perturbations at the relatively low redshifts at which type la supernovae have been observed. 

An initially Gaussian spectrum of density perturbation will evolve into a log-normal 
probability distribution of the form ^Hl 

PiS)dS ^ exp (ZM^±3±^) ,s . (18) 

where the density contrast 6 = Sp/p. The parameter a has been calculated numerically 
using simulations of cosmologies with and without a cosmological constant ^T^. We are not 
aware of such a calculation for a string network cosmology but since the variation between 
the sigma for a flat universe with or an open universe without a cosmological constant is of 
order 20% [TT], we proceed using the expression for the = 0.3, I^a = universe: 

a{z) = 0.08096 + 5.3869(1 + z)-^ - 4.21123(1 + z)-^ + 1.4433(1 + z)-\ (19) 

One should keep in mind that with this density distribution, most cells are under-dense 
relative to the mean by a factor of about ~ 10. 

Since the number density of the electrons traces that of matter, and the average num- 
ber density of electrons goes as the third power of the size of the universe, we obtain the 
relationship 

n,{z) = il + 6iz))il + zffi,{0). (20) 

This quantity is calculated in each cell, taken to have the size of the Jeans length, along the 
line of sight using the the probability distribution Eq. ()18|1 . For the magnetic field, being 
frozen into the medium, we have the relation 

B^{z) = (1 + 6{z))^ (1 + zf B{0) sinV' (21) 

where is the angle between the magnetic field and the line of sight. Since there is more solid 
angle perpendicular to the line of sight than along it, we took into account the appropriate 
trigonometric weighting when randomly generating ip for each cell. Many uncertainties 
persist regarding the structure of intergalactic magnetic fields and their origin is a 
much debated topic PH]|in]- Here we simply assume that the magnetic fields were created 
before any of the redshifts that we are interested in. Also, we approximate the coherence 
scale of the magnetic field to be of the order of the Jeans length. 

Now we have set up the background cosmology, defined the equations for mixing as a 
function of electron density and magnetic field and determined how these two components 
behave with redshift, we proceed with the simulations. 



3 Results 

We performed simulations to calculate the effect of axion dimming upon light from high 
redshift objects based upon the cosmological parameters discussed in the previous section. 
We also investigated the energy dependence of this dimming and the effect of the mixing 
upon the polarisation of distant galaxies. 
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3.1 Photon Dimming 

If we consider some flux of light travelling through a cell with a certain probability of mixing, 
we need to see what proportion of that flux will convert into axions at the end of its path 
between source and detector. As mentioned earlier, since the coherence of the magnetic field 
is lost from cell to cell have taken the orientation of the magnetic field in each cell to be 
random, implying that axions will on average mix into photons with the same probability. If 
we use / to denote the fraction of the original flux of light from some source which remains 
in the form of photons we can find the value of / after some time tn in terms of the mixing 
matrix p(tn) 

/(tn) = (22) 
^j=l,2,3Pjj[tn) 

which comes about because there are two polarisations of the photon but only one axion. 
We simulated random lines of sight out to redshifts of 2; = 0.5 and z = 1 since these are the 
redshifts around which the type la supernova observations are concentrated [HEl- 

The mass-density within each cell was randomly generated with weighting from the den- 
sity probability distribution (fTHj). then equations PTj) and were used to obtain the 
strength of the magnetic field and the electron density. The size of each cell was set to the 
jeans length corresponding to the density of that cell. We set the overall magnetic field 
strength and electron density by defining its average value at 2; = 0. The electron density at 
2; = is fixed to be ne(2; = 0) = 1.8 x lO^'^cm^^ by observations of the primordial abundance 
of light elements whereas we are have more freedom in choosing our value for the magnetic 
field. Since the effect of changing the photon axion coupling Ma and the magnetic field leads 
to a degeneracy in our results, we instead choose to use define a dimming parameter t] as 

/ B(z = 0)\ flO^^GeV\ 1 , , 



The colour excess E{B — \^) is a parameter normally used to quantify the difference in 
the ratio of emitted and observed flux (denoted Femi and Fobs respectively) between the B 
and V bands due to intervening matter. It therefore gives an indication of the frequency 
dependence of any dimming between the observer and the source. The colour excess is 
defined as 

Supernova observations have placed an upper limit upon this quantity of E{B — V) < 0.03 
[5j so we must try to vary 77 to obtain enough dimming in order to explain the supernova 
results whilst staying below this colour bound. Figure [T] plots the quantity E{B — V) against 
the dimming parameter rj based upon our simulations using the background cosmology and 
density probability function described earlier. This figure shows that the maximum amount 
of permitted dimming at a redshift of z = 1 which would still satisfy the colour bound is 
about 12 %. In figure |21 we have plotted photon dimming vs. the parameter 77. A dimming of 
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Figure 1: Colour excess vs. dimming of photon flux due to mixing into axions for different 
values of rj at redshifts z = 0.5 and z = 1 
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Figure 2: Dimming of photon flux vs. rj at redshifts z = 0.5 and z 



12 % at 2; = 1 indicates a value for this parameter of ~ 1 x lO^^^cm^^ which also suggests 
a dimming of only about 8% at 2; = 0.5. A larger dimming which would better account for 
the data at 2; ~ 0.5 is rj = 1.5 x 10~^^cm"^ although one might then have already expected 
to see supernovae which violate the colour bound. 

3.2 Polarisation 

Similarly we can investigate the evolution of light from a polarised source (see also [20112] for 
related previous work). The magnetic flelds in between us and the polarised source will not 
of course be aligned with the polarisation axis of the source. Photons will therefore be able 
to mix into axions and axions into photons and the net angular momentum of the flux will 
not be conserved since the mixing between the photons and the axion fleld transfers angular 
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momentum to and from the surrounding magnetic field. By carefully analysing the mixing 
we are therefore able to find out the way a highly polarised source loses its polarisation as 
the photons propagate through the intergalactic medium. As in the case of dimming, it is 
important for the magnetic field orientation and the direction of propagation of the photons 
to be uncorrelated in order for information concerning the polarisation to be lost in this way. 
However, since we expect there to be a large number of cells with randomly oriented B-fields 
in between us and the source, this does in fact occur. 

Let us consider a completely (100%) linearly polarised source at a particular redshift. 
As before, we call the two axes perpendicular to the direction of photon propagation and 
each other x and y. We choose to align x with the polarisation axis of the source which 
corresponds to an initial density matrix of the form 



and then calculate the mixing of the photons from such a source. As these photons pass 
through the inhomogeneous matter distribution in our simulations we can define a polarisa- 
tion ratio n at a time t„ after n cells as 



which tells us how much the linear polarisation of the light will have been diminished in 
reaching us. 

Figure 121 shows the observed polarisation as seen on earth of a completely polarised source 
located at a redshift z assuming the light from the source has been diluted by photon-axion 
mixing along the way with a dimming parameter of rj = 1.5 x 10^^^ cm^^. We have considered 
mixing of light from four different parts of the spectrum, the B and U bands in the optical, 
the boundary between the 'visible' and the ultraviolet window (300nm) and well into the 
ultraviolet (lOOnm). Any astronomically observed source with redshift and polarisation 
values which lie above these curves will rule out this degree of mixing. 

Tabled contains a few results of polarimetric observations of high redshift active galactic 
nuclei^. The values listed in this table are compatible with this value of t], so it would not 
be necessary to reduce the amount of dimming in order to accommodate these observations. 
Consequently, to our knowledge, polarimetric observations do not yet rule out cases where 
the dimness of supernovae can be explained by axions. 

However, the polarisation effect is extremely frequency dependent, and figure El shows 
that we do not expect to see any highly polarised high redshift sources in the ultra-violet 
part of the spectrum (A ~100nm) with values of rj which could lead to significant dimming. 
UV Polarimetry of high redshift objects must be done from space and table ^ includes 
observations using the Hubble space telescope faint object spectrophotometer [2^1 ■ Such 
observations should also be possible with the new instruments on the Hubble space telescope. 

^ Tabled does not represent an exhaustive literature search so it is entirely possible that data exists which 
could constrain the parameters more strictly. 




(25) 




(26) 
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Figure 3: Observed polarisation ratio 11 of completely polarised sources located at different 
redshifts z and seen at different A. Solid line is B band (440 nm), dashed line is U band (360 
nm), dotted line is edge of UV spectrum (300 nm) and dot-dash line is 100 nm. Dimming 
parameter 77 = 1.5 x 10~^^cm~^. 



Table 1: Some observations of highly polarised high redshift active galactic nuclei. \ohs is 
the central wavelength of the light used to obtain 11 in the rest frame of the telescope. 



Object 


Redshift (z) 


Ao6^(nm) 


n(%) 


Reference 


4C+03.24 


3.56 


604 


11.3+3.9 




22 




4C+00.54 


2.37 


443 


11.9+2.6 




23 




0211-122 


2.34 


435 


19.3+1.1 




22^ 




FSC 10214+4724 


2.28 


416 


26 +2 




24 




0823-223 


>0.91 


350 


16.8+0.84 




25^ 




1522+101 


1.32 


166 


4.7+1.9 




26, 
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4 Discussion 



In this paper we have considered the dimming of hght from high redshift supernovae due 
to mixing between photons and a pseudoscalar axion field in the intervening intergalactic 
plasma. We have tried to simulate this process in a more detailed way than previous studies 
by taking a realistic density probability distribution for the electron density. We have also 
included the effect of a fluctuating magnetic field due to it being frozen into the background 
plasma. 

The presence of such an axion does not completely solve the problem of the missing ~ 70% 
of the energy of the universe. Rather we aimed to show that such a mixing might make it 
possible to explain the missing dark energy of the universe via an alternative source of stress 
energy other than a cosmological constant. As a candidate for this energy we considered the 
case of a string network with equation of state us = —1/3 and calculated the combination 
of axion-photon coupling and magnetic field that would lead to enough dimming. 

Assuming a primordial origin for the intergalactic magnetic fields we found no values 
of these parameters which could explain the dimming without violating the bound on the 
frequency dependence of this dimming already calculated by the observational supernovae 
groups. However, the colour bound that we observed was only a factor of about 50% larger 
than that observed, not many orders of magnitude, and there is a small inherent error 
expected in our analysis due to deviations from our assumed behaviour of the variance 
parameter a{z) in the density probability function ()18|) which will become greater at high 
redshifts. We found that for rj = 1.5 x 10~^^cm~^ we obtained nearly enough dimming and 
only marginally violated he colour bound. We therefore have adopted this as our optimal 
value of rj. 

One of the main conclusions of this work is therefore that a photon-axion coupling of 
the sort discussed in [HI leads to very particular predictions as to the redshift dependence of 
dimming at different frequencies. 

If one considers an alternative origin for the intergalactic magnetic fields other than their 
primordial production it will be possible to change the bounds obtained in this paper since 
the correlation between magnetic field strength and electron density used in our calculations 
will be modified. 

We also calculated the expected reduction in the polarisation of light from high redshift 
objects due to photon-axion mixing. After a brief literature search we were unable to find 
any observations of high redshift AGN with large enough polarisations to be at odds with 
a dimming of 77 = 1.5 x 10~^^cm~^. However, we find the frequency dependence of this 
effect to be so large that we believe such a dimming would lead to no observed sources 
at redshift z ~ 1.5 with UV (100 nm) polarisation as high as 10%. Any distribution of 
magnetic fields which could lead to significant dimming of the optical flux coming from a 
distant object should also lead to a large reduction in any polarisation of the UV photons 
coming from that object. This prediction is independent of the detailed history of cosmic 
magnetic fields. It should be possible for this issue to be further investigated by making 
space based polarimetric observations. 
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